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SUMMARY 



The first step in the development of a hybrid com- 
puter simulation of the data-based, 5 DOF, 3K-SES program 
that will operate in real-time has been completed. A 
3 DOF (Flat) turn model for an approximate and simplified 
XR-3 craft was used to provide the equations of motion in 
order to become familiar with 1) the turning character- 
istics of the CABSES, 2) the requirements of the (NPS) NAV 
PG School hybrid computer with regard to signal conversion, 
timing, accuracy and loop delay due to both signal con- 
version and (right-hand-side) computation and 3) the 
approximate size of program that can be computed on the 
existing hardware of the NPS facility. 

In parallel with this investigation, the 5 DOF data- 
based digital program for the 3K-SES has been converted 
from CDC to IBM FORTRAN and is now operational and will 
serve as the standard for evaluation of the final hybrid 
simulation for the 3K-SES . 

Results of the hybrid computer simulation of the 3 DOF 
turn model of the XR-3 indicate that the installation of the 
5 DOF program of the 3K-SES should proceed (in the next 
phase) as a straight forward application of the right-hand- 
side computations within the digital computer and addi- 
tional degrees of freedom on the analog computer. The 
successful operation of the hybrid computer simulation will 
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depend upon ability of the 10 year-old hardware--in 
particular the digital computer--to compute the many 
computational operations required for the right-hand- 
side of the equations of motion in the time frame necessary 
to approximate real-time manual control. 
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INTRODUCTION 



* 
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Reference 1 presents the Rohr Marine, Inc. analysis 
of the stability and maneuverability of the 3K-SES under 
direct manual control. The recommendations contained in 
Reference 1 included a suggestion to continue the investi- 
gation of manual control using a more realistic model of 
the response of the human operator to emergency conditions. 

If a real time simulation can be developed it can 
provide a means of inserting an actual human operator in 
the loop, as well as any selected model of the operator. 

A hybrid realization of such a real time simulator appears 
to be feasible and would provide a convenient hands-on 
facility for a variety of human operator studies, plus 
many other studies in control and guidance. 

This report presents the results of the first phase 
of a feasibility study, in which it is planned to develop 
a real time simulation using the NPS Hybrid Computer 
which consists of a CI-5000 analog unit and an XDS-9300 
digital unit. The initial studies reported here are 
concerned with development of 3 DOF models for the XR-3 
and the 3K, as a preliminary to the development of a 5 DOF 
Hybrid model. 
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II. FAMILIARIZATION AND PLANNING 
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IIA. THE 3 DEGREE OF FREEDOM MODEL 



1. INTRODUCTION 

The development of any computer model requires famili- 
arity with both the system equations and the computer to be 
used. In addition, one needs quantitative data to insert 
the model into the computer, and one also needs typical 
performance data such as experimental test results in order 
to validate and verify the model. To develop a hybrid model 
of the 3K-SES, we chose to start by modelling the XR-3, 
because we have the required quantitative data, we have a 
good 6 DOF digital program, and we have good experimental 
data, plus availablility of the XR-3 craft for test as needed. 
We also decided to start with a 3 DOF model, for several 
reasons : 

a) It is the simplest model which will permit turning 
tests that are verifiable. 

b) The complexity of the initial model will thus be 
minimal, but we can readily add two more degrees 
of freedom with confidence in the resulting 5 DOF 
model, because the similarity of the differential 
equations will permit ready application of the 
techniques we develop in perfecting the 3 DOF model. 

In addition to this decision to start with a 3 DOF model, 
we chose to simplify the model by simplifying the right-hand- 
side of the equations . Since the terms on the right-hand- 
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side are arithmetic calculations to be done in the digital 
unit, and because the many terms are additive , truncation 
of the right-hand-side, with simplification of the terms 
retained, should have little effect on results except some 
decrease in quantitative accuracy. This permits us to work 
with a less complex model in the initial stages, and addition 
of terms at a later point should be a simple matter. It 
will also have the advantage of permitting us to assess the 
effects of adding terms on the computational delays and on 
the limits of capability of our hybrid installation. 

2. EQUATIONS OF MOTION OF THE 3 DOF MODEL 

The three degrees of freedom to be used in this devel- 
opment are: SURGE, SWAY, YAW. As a result the following 

assumptions and approximations are needed: 

1. Roll is zero 

2. Pitch and draft remain constant 

3. The origin of the ship’s coordinate system is 
located at the center of gravity. 

4 . At t = 0 the ship coordinate system and the fixed 
(geographical) coordinate system have the same 
origin and the same coordinate axes. 

5. At t = 0 the ship is heading in the positive 

X „ direction, 
nav 
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Figure 1. DEFINITION OF COORDINATE SYSTEMS. 
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The equations of motion can then be approximated by: 



SURGE m(u - vr) = X (1) 

SWAY m(v + ur) = Y (2) 

YAW I r = N (3) 

Navigation 

x Q = u cos V - v sin ¥ (4) 

y Q = u sin y + v cos V (5) 

v = “\|u+"”v^ (6) 

5 

3 = tan" 1 (7) 

In the above equations: 



m = mass of the rigid ship 

I = moment of inertia about the z-axis 
z 

u = v g cos 3 = velocity in the x-direction (SURGE) 

v = V sin 3 = velocity in the y-direction (SWAY) 
s 

r = ¥ = angular velocity about the z-axis 
y = heading (yaw) angle (see Fig. 1) 

3 = drift angle (see Fig. 1) 

5 = rudder angle, + <$ produces a left turn (see Fig. 1) 



Also: 

x o = A dt 

i’o = / : 'o dt 
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The right-hand-sides (REIS) of the equations of motion 
have been designated X, Y, N. These symbols represent 
summations of forces (X,Y) and moments (N) , and each summation 
may consist of many terms. For the initial simplified model 
we have discarded all terms in the summations except the 
few that are obviously major contributions. The terms re- 
tained are considered sufficient to produce a working 
simulation which should provide qualitative agreement with 
reference data. At a later point the model may be refined 
by adding terms and adjusting terms and adjusting coefficients. 



After reduction of the right-hand-sides, these force 
summations become: 



X 


2 

= GT cos 6 - C v u 


(8) 


Y 


= GT sin 6 - C_ vlvl 
Dy 1 1 


(9) 


N 


= -(GT sin 6)1 + vlvll 

o Dy 1 1 a) 


(10) 


where : 


(see Fig. 2 also) 





GT is the gross thrust in pounds 

<5 = <5 =6 = angle of all thrust vectors 

s P 

C_ = coefficient of drag in the direction of surge 
D-K 

= coefficient of drag in the direction sway 
1 (see Fig. 2)= thrust lever arm 
1 (see Fig. 2)= sway-drag lever arm 

In further explanation of equations 8, 9, 10, note 
that we have chosen to model the drag forces as lumped 
parameters acting at two points P and P ” (see Fig. 2) . 
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2 

NOTE « 8 = 8„ = 8 

$ p 

P - P 1 ARE EQUIVALENT FORCE CENTROIDS 
Figure 2 . FORCES ACTING FOR A RIGHT TURN. 
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The length 1 is the distance from the stern to the center 
of gravity, and 1 is the component distance (in the x- 
direction) from the center of gravity to the points P-P'. 

1 is a fixed parameter for our initial studies (it could be 
changed by ship loading) but 1 is an adjustable distance 
which we will use as a parameter to make the model fit test 
data. 

The coefficient C DX is a lumped parameter equivalent 

for the skin friction effect (slender body theory) . The term 
2 

C D ^u models the hull drag, aerodynamic effects and hydro- 
dynamic effects. The bubble drag is not included in the 
initial model (will be addded later) because our initial 
efforts are aimed at evaluation of our hybrid structure, 
and a simple model above critical speed will suffice for 
this. In like manner, in eqn. 9, the coefficient C D y models 
the cross flow drag effects. Note that the added mass effect 
has been temporarily omitted since its effect is also of 
second order. 

Substituting equations 8, 9, 10 into equations 1, 2, 3 
provides the basic equations of motion for the 3 degree of 
freedom model: 

u = (GT cos 5 - C^u 2 ) + vr (11) 

v = (GT sin 6 - C Dy v|v|) - ur (12) 

r = (" (GT Sin 6)1 o + (C Dy v|v l )1 a J ) (13) 

Equations 11, 12, 13 are applicable to any SES craft. 
Selection of parameter values to be used in the equations 
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would depend on the particular craft, but the procedures 
used to identify the lumped parameter values are independent 
of the craft type. 

3. IDENTIFICATION OF CRAFT PARAMETERS 

In order to evaluate the parameters needed for the equa- 
tions, the craft must be available for measurement and test 

(or suitable data must be available) . The parameters m, I , 

z 

1 are constants which must be known, and GT, 3, u, v, r 
can be measured in steady state (if needed) . Then C DX , 1 , 
C D y can be calculated as is shown below: 

A. Surge Drag Coefficient, C DX 

The craft is run on a straight course with GT = 
constant and given rudder conditon 5 . When steady state is 
reached u is measured. Then 

i = 0 = (GI as J - (14) 

fram which 

C DX = GT COS 6/u ss U5) 

where 

GT = gross thrust in pounds 

u ss = steady state speed in feet/sec 

5 = rudder angle 

B. Sway Drag Moment Arm, 1 

Let GT = constant and 5 = 5^, a fixed angle. When 
the craft reaches steady state in the turn the sway equation 
becomes 

v = ° = — (GT sin 3 - C v|v|)- ur (16) 
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from which 



C_ v 
Dy 



= GT sin 5 - m u r 



The steady state yaw moment equation becomes 



r = 0 = 



[- 



(GT sin m 0 + (C Dy v|v|)l u 



from which 



] 



C v 2 = (GT sin 6) ( 1 o/ 1 ) 
^ co 

Combining equations 17 and 19: 



(17) 



(18) 



(19) 



GT sin 5(1 - °/ l ) 



r = 



mu 



or if we wish 
1=1 



GT sin 5 



( 20 ) 



( 21 ) 



o GT sin 5 - m u r 
From Fig. 2 and eqn . 20 it is clear that the turning 
rate, r, is a function of the distance 1 , with all other 

OJ 

parameter values fixed. Since it is customary to command 
a turning rate, if we choose r we can then evaluate 1^. For 
the XR-3, using steady state data obtained with our 6 DOF' 
simulation , 

r = -0.0266 (1 - ) radians/sec (22) 

co 

A plot of this relationsip is shown on Fig. 3. 

C. Sway Drag Coefficient, C D 



Again assuming steady state conditions with fixed 

thrust and fixed rudder, the sway equation, eqn. 17, provides 
^ -m u r + GT sin 5 



'Dy 



(23) 



v 
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Figure 4-A. FLAT TURN DYNAMICS OF THE XR-3, PLOT 
OF VELOCITY vs. TIME, 
r = .33 RAD/SEC ( = 1 ft) 

W 
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Figure 4-B FLAT TURN DYNAMICS OF THE XR-3, PLOT 
OF X vs. Y. 

r = .33 RAD/SEC (0 = I ft) 

w 
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Figure 5-B FLAT TURN DYNAMICS OF THE XR-3 PLOT 
OF X vs. Y. 

r = 0.1 RAD/SEC ( ft =2.5 ft) 

w 
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For a given turning rate, r, we know all values except v, 
which can be calculated using the surge equation, eqn. 11, 
to provide: 



v 



C^u - GT cos 5 
_ m r 



(24) 



4. DIGITAL SIMULATION OF THE SIMPLIFIED XR-3 EQUATIONS IN 
3 DEGREES OF FREEDOM 



The equations of motions as previously derived were 
programmed into the IBM 360/67 computer. A fixed step size 
was used to verify the program and validate the model. 

Using a turn rate, r = .33 rad/sec, and an integration step 
of 0.1 seconds the simulation was run for 120 seconds of 
problem time, then repeated for an integration step of 0.01 
seconds. Comparison of results showed that an integration 
step of 0.1 seconds was suitable, all variables in the two 
runs agreeing to within 1% . 

Figure 4A and 4B show plots of u vs time and x vs y 
for a flat turn with r = .33 rad/sec (1 = 1 ft.) . Figures 

U) 

5A and 5B show similar plots for r = 0.1 rad/sec (1 = 2.5 ft.). 

A number of additional runs were made but results are not 
given here. Based on these runs we feel that the program 
has been satisfactorily verified, i.e., it is solving the 
equations correctly. 

The following observations seem appropriate: 

1. Eqn. 22 and Fig. 3 are at least qualitatively 

correct . 

2. The effect of added mass is neglegible. 
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3. We have gained confidence that, given the meas- 
ured performance of any SES in an approximate flat turn, and 
given values for the parameters GT, 5, u, r, m and 1 q , we 
can determine a value of 1 which provides a simulation 
matching the measured performance reasonably closely. 

4. If the simulation is performed with an initial 
step in the thrust angle, <5, approximate values for the system 
time constants can be evaluated from the resulting data. 

For the XR-3 operating in a flat turn the time constants are: 

Variable Time Constant 

u 11.0 sec 

v 0.8 sec 

r 1.2 sec 

0 0.9 sec 

To validate the 3 DOF model, additional runs were made 
using the digital 6 DOF Oceanics program as modified by our 
NPS group for the XR-3. Also the 3 DOF model was run with 
3K-SES data, and the 5 DOF model of the 3K-SES was run. For 
the XR-3 tests a step of rudder of 0.5 radians was suddenly 
applied with craft speed at 20 knots. For the 3K-SES tests 
and 11° rudder was applied by ramping up for 1.0 second and 
holding at the 11° value. Craft speed was 60 knots. Com- 
arison of the 3 DOF, 6 DOF and 3 DOF Hybrid runs for the 
XR-3 is shown on Fig. 6. Correlation is not as bad as 
appears by inspection. Turning diameter and time for a 360° 
turn are almost indentical. We presently attribute the 
differences to the fact that the 3 DOF model does not include 
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Figure 6. TURNING CHARACTERISTICS OF THE XR-3 RESPONSE 
TO LEFT RUDDER, 0.5 RADIAN STEP, 20 KNOTS. 
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a term for bubble drag. 

Figure 7 compares the 3 DOF run with the 5 DOF run 
for the 3K-SES. The results are almost identical, despite 
simplifications in the 3 DOF. This gives confidence that 
our modelling technique will work well when we hybridize. 

5. HYBRID SIMULATION OF THE 3 DEGREE OF FREEDOM MODEL OF 

THE XR- 3 

In developing a hybrid model various arrangements and 
structures (programming) are possible. The philosophy 
chosen for our initial hybrid model is shown in the block- 
flow diagram of Fig. 8. Details of the programming are 
given in Appendix A for the simplified XR-3 model. 

The simplified XR-3 model was installed in the hybrid 
system and tested. Comparison of results with those ob- 
tained by digital simulation showed almost identical per- 
formance, differences being less than 2%. 

In developing the hybrid simulation a major concern 
has been the amount of time delay introduced into the loop 
by the A/D and D/A conversions and by the time required to 
perform the digital computations. We have monitored these 
delays in the 3 DOF model with results as follows: 

1. The fundamental loop delay using D/A for u, v, r, 
x, y and A/D for u, v, r was 33 ms. 

2. A lower language subroutine exists, called ADDA, 
which can be used to replace the Fortran subroutines 
for A/D and D/A. This can reduce the delay by 25%. 
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Figure 8. STRUCTURE FOR HYBRID SIMULATION OF 3 D. 0. F. (FLAT) TURN 
EQUATIONS. 







3. Additional delays are introduced by using A/D 

for thrust (GT) and thrust angle (5) C 76 % increase] 

and by use of the Graphics Signal (Gated) conversion 
of the variables [6 7% increase] . Total loop trans- 
port delay using D/A, both A/D units and the Graphics 
signal conversion was approximately 7 8 ms . 

We know how to reduce the loop delay substantially, but 
in view of the system response time this may not be necessary. 

One possible source of difficulty in the hybrid simulation 
is the cross coupling between the equations . When operating 
near steady state conditions during small yaw angle maneuvers 
the linearized equations indicate that the roots of the 
characteristic equation may be very lightly damped. 

Consider straight-ahead motion (thrust) at constant 
speed. The linear model sway and yaw equations are: 



*2 = 


a 22 X 2 " U o X 3 


*3 = 


a 32 X 2 


'' A 

X = 


|jAu, iv, Ar"j T , and Au = 0 


a 22 


f l (x 2 ) 


a 32 = 


f 2 ( x 2 ^ 



Combining these equations, the characteristic equation 

of the simplified 3 DOF model is 
2 

s + a„ 0 s + u a,- = 0 
22 O 3 2 

The locations of the roots of this equation depend on the 
values of a £2 and a 22' koth °f which are functions of the 
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sway velocity. For small values of sway velocity the 
roots are complex conjugates and approach the origin of 
the s-plane as the sway velocity approaches zero. For 
larger values of sway velocity such as those encountered 
in a turn maneuver the roots are real and negative, pro- 
viding a well damped response for yaw and sway motions. 

This appears to be similar to the Dutch Roll effect in 
aircraft. 

For straight ahead runs using the digital program for 
the 3 DOF model and setting the initial value of v to zero, 
no sway-yaw motions are observed. However, if we set a 
small initial condition on v a transient oscillation is 
observed in both sway and yaw as shown on Fig. 9. This 
verifies the existence of the phenomenon, and we expect 
to encounter it in the hybrid since the equations are the 
same. However, in the hybrid there is a possibility that 
the roots may move into the right half s-plane due to either 

1. DC unbalance in the analog computer integrators, or, 

2. Excessive loop delay from A/D, D/A and digital 
computations . 

The first of these possible causes can be avoided by 
careful balancing of the operational amplifiers before 
running. The effect of the loop delays can be controlled 
somewhat by minimizing the number of terms requiring digital 
comput at ion. 

6. STATUS OF THE 3 DOF MODEL STUDIES 

These appear to be both complete and successful. Both 
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qualitative and quantitative agreement have been achieved 
with the 5 DOF digital simulation as a reference. Future 
studies will be aimed primarily at the 5 DOF model of the 
3K-SES except when we wish to correlate computer data with 
specially designed experiments. 
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IIB. THE 5 DEGREE OF FREEDOM MODEL 



1. INTRODUCTION 

Having developed a simplified 3 DOF model of the XR-3, 
we now proceed to use the knowledge thus obtained and apply 
it to the 3K-SES. We do not plan to develop a 5 DOF 
hybrid model of the XR-3. Our procedure is to install 
the NSRDC 5 DOF program for the 3K-SES in our IBM 360/67 
computer, generate reference curves with it, obtain any 
numbers we need for our 5 DOF hybrid simulation, then 
verify and validate the hybrid model of the 3K-SES. 

2. THE DIGITAL SIMULATION OF THE 3K-SES IN 5 DEGREES OF 

FREEDOM 

This program was obtained from NSRDC. It has been 
written for a CDC machine and had to be converted to the 
IBM version of Fortran for our computer. This has been 
completed, the program is operational and as far as we 
know has no obvious "bugs”. We are familiarizing ourselves 
with the program, and exercising it to obtain data we will 
need in developing the 5 DOF hybrid model, and for validating 
it. 

Details of the conversion-problems encountered, their 
solutions, etc. are discussed in Appendix B. 

3. HYBRID SIMULATION OF THE 3K-SES 

This is under development at the present time. Hope- 
fully we need only insert a new set of constants in the 
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equations of the hybrid XR-3 model, but undoubtedly other 
adjustments and modifications will be needed, particularily 
in the trade-off of digital vs analog computation in order 
to perform a real-time simulation. 

4 . STATUS 

Progress has been satisf actory . It is too early to 
predict future developments. 
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III. SUMMARY AND CONCLUSIONS 
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The NPS Hybrid Computer consisting of the XDS-9300 
digital and the Cl 5000 analog computers has been pro- 
grammed to simulate the 3 DOF Flat Turn maneuvering 
dynamics of the CAB type Surface Effect Ship. 

The XR-3 has been programmed with simplified modelling 
of the drag forces and thrust-input force on the hybrid 
computer for simple flat turn maneuvers. An evaluation 
of the hybrid simulation was carried out by comparison of 
the time histories of the velocities and the navigation 
variables against those values obtained from the 3 DOF 
digital simulation for the same maneuvers. In addition, 
the XR-3 hybrid 3 DOF results were compared to the 6 DOF 
loads and motion output . 

The result of this comparison was that the hybrid 
computer simulation produced a time response to a step 
rudder input that produced a navigational turning plot that 
was in-between the all digital 3 DOF and all digital 6 DOF 
trajectories . 

Additionally the 3 DOF simplified RHS (right-hand-side) 
digital program was operated for the same turning conditions 
that was used for the 5 DOF, DRSIM5D Program and the 
comparison indicates close agreement in the navigational 
trajectories for an 11 degree vectored thrust angle full- 
turn maneuver. 
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This feasibility study has shown that the hybrid 
computer implementation for the 5 DOF 3K-SES in real- 
time is possible with the available equipment at the NPS 
Computer Laboratory. Real-time computation is of course 
only approximated when using the hybrid computer. The 
loop transport delay inherent in this simulation is due 
to the signal conversion (A to D and D to A) and the 
necessary digital computation for the acceleration forces 
and moments. The magnitude of this delay time has a fixed 
minimum value as a result of signal conversion and a maximum 
value due to righ-hand-side calculation required for the 
acceleration added to the fixed time due to the signal con- 
version. Real-time computation accuracy is determined 
principally by the magnitude of this maximum value for the 
loop delay. 

The 3K-SES was approximated with the simplified 3 DOF 
Flat Turn equations and the results were compared. It 
was shown that the full maneuvering trajectories for an 
11 degree thrust angle compared very closely with the 5 DOF, 
DBS IMS D Program results. These results indicate that only 
a few terms on the RHS may be all that are required to simu- 
late the 3K-SES for the type of maneuvering control studies 
to be conducted in the future on the hybrid computer. 

The conclusion of this feasibility study is that suc- 
cessful implementation of a hybrid computer, real-time 
simulation of the DBSIM5D Program has a very high prob- 
ability. The degree of complexity necessary for the RHS 
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computation will be the limiting factor on the loop delay 
that controls the real-time accuracy. 
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APPENDIX A. HYBRID COMPUTER PROGRAM FOR THE 



3 DOF (FLAT) TURN EQUATIONS 

A general schematic diagram of the hybrid structure 
used is shown in Fig- A-l. The listing of the digital 
program for the XDS-9300 computer is also included in this 
appendix, and the analog computer connections are shown 
on Fig. A-2. 

As may be seen from the listing of the digital pro- 
gram, the right-hand-sides (RHS) of the acceleration 
equations for the 3 DOF are computed on the XDS-9300, 
which also does the magnitude scaling for the analog 
computer (CI-500) voltages. In addition the XDS-9300 
provides the required control logic, timing and temporary 
storage of computed variables. 

The graphics terminal (AGT-10) provides real time 
display of the desired output and control variables, as well 
as their maximum and minimum values. Piloting of the craft 
is carried out at the graphics terminal. Additional 
permanent record, outputs of navigation x and y are obtained 
with an XY-pen recorder, and plots of other system variables 
are obtained on an 8 channel Brush recorder. 

A sample of the x-y output for the 3 DOF simplified 
XR-3 craft is shown on Fig. A-3. The plot starts with a 
right turn (6 = -0.5 radians) at t = 0 , using a terminated 
ramp input. The rudder was held at 6 = -0.5 for two turns, 
then reversed to 6 = +0.5 for a half turn, followed by 
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HYBRID COMPUTER 




Figure A-1. BLOCK DIAGRAM OF THE HYBRID STRUCTURE. 
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Figure A- 2. DIAGRAM OF THE ANALOG COMPUTER 
CONNECTIONS. 
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<5 = 0.0, then again 5 = +0.5 for a quarter turn, and finally 
a right turn using 6 = -.29 radian for two complete turns. 

The total time for these maneuvers with the XR-3 was approxi- 
mately 3^- minutes. 

Note that the CI-500Q is used only for analog inte- 
gration of the accelerations to obtain the desired velocities 
for the digital computation. Displacement is computed in 
the digital computer using the trapezoidal rule for integratior 
as done on the Oceanic's 6 DOF Loads and Motions Program. 

Note also that the present structure and programming 
use only existing hardware and we consider them a feasi- 
bility study. If results show that we have been successful, 
then we must prepare for manual control studies, and this 
may require development of a more realistic ''cockpit” to 
simulate conditions in the control room of the 3K-SES. 
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Figure A-3. SAMPLE OF THE HYBRID 
OUTPUT x-y PLOT FOR 
THE 3 D.O.F. XR-3 
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APPENDIX B. INSTALLATION OF THE 3K-SES , 5 DOF 



PROGRAM IN THE IBM 360/67 COMPUTER 

In the fall of 1978 the 3K-SES , 5 DOF data based 
program called DBSIM5D was delivered to the Naval Post- 
graduate School, at which time a commitment was made by 
the NPS group to attempt a hybrid computer simulation based 
on this program. 

In the spring of 1979, DTNSRDC shipped the card deck 
and the output data for a typical 3K-SES run, obtained by 
DTNSRDC on a CDC computer. The deck was to be used in 
installing the program in the IBM 360/67, and the output 
data provided information for verification and validation. 
The major problem in the process of installation was the 
conversion from CDC Fortran to IBM Fortran. With the help 
of Mr. Roger Killeary of the NPS Computer Facility, the 
interactive graphics terminal was used to edit the program. 
After some 50 man hours of effort the first successful 
batch run was made in April 1979. 

Complete success in the conversion did not occur until 
a correction was added to the newly programmed TIMER, to 
allow for a stop at the desired stop time rather than a 
premature stop time due to the timer program. 

To validate the IBM version of the program we dupli- 
cated the sample run for the 3K-SES . This run, as provided 
by DTNSRDC, was for a left turn, using an 11 degree effector 
angle and a speed of 60 knots. The sample run terminated 
after 14.5 seconds of real (problem) time. The 11 degree 
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effector angle was programmed as a ramp which terminated 
after 1 second, and remained constant at 11 degrees for 
an additional 50 seconds. 

A comparison of the output of the IBM 360/67 version 
with the sample run data from the CDC computer showed that 
the recorded variables were in complete agreement within 
the accuracy differences of the two computer word length. 

The problems of conversion were simple but time 
consuming and may be listed as follows: 

1. 6 character variable name limit on the IBM com- 
puter as compared with, a 7 character limit on the 
CDC. 

2. Single statement limit on the IBM computer versus 
multiple statement on a single card for the CDC. 

3. Identification of temporary files used for the 
output data. 

4. Identification of the timing control and pro- 
gramming the IBM computer timer to provide the 
functions required for operation of the program. 

5. Addition of the proper job control language (JCL) 
cards to compile, link and compute. 



DAT 



Are D TR A \ 

l : 
o • 

3 : 
a : 

5 ; 
6 : 

7 ; 

5 : 

Q ‘ 

io: 

1 1 : 
13: 
13: 

1 u : 
15: 
16: 
17: 
13: 
19: 
30: 
21 : 
32: 
23: 
24 : 
25: 
26 : 
27 : 
22 : 
29 : 
30 : 
31 : 
32: 
33: 
34; 
33 : 
36: 

37 ; 

38 : 
39: 
40 : 

4 1 : 
42 : 
43 : 

44 : 

4 5 : 
<*6 : 
47: 
4? : 

49 : 
50 : 

5 1 : 



LS, GP 

Cl YENS 1 9N ITrXT(29),ILAP(24),P9 r 3(24),!=5TK(24), I * AGE ( 3003 ) 

C, I L ( 24 ) , TITLE ( 24 ) , AC < 10) , AA ( 12 ) , ILA1 (24 ) , ILA2(24) , I „a3( 24) , 
CILA4(04)/ ILA5(24)^ ItA6(24) / IL*7(24) , ILAS(24) / IRt.3T< 10), IvIEaA(2 
C I V I EwE ( 2C ) 

EQUIVALENCE (AD(1)*ai1),(AD(2)#A21).(AD(3)/A31)/(AA(1)j311)* 

C< AA(2) ,321 ) , ( AA( 3), S31 ) , ( AA( 4 ) ,341 ) , ( AA (5) ,351 ) » ( AA< 6), 355 ) , 

C( AA(7 ) ,361 ) , < AA( 8) , 37 1 ) i ( A A (9 ) , 331 ) , ( A A ( 10) ,391 ) , ( A A ( 11 ) , 3101 ) , 
C( AA < 12 ) ,B11 1 ) , ( A0( 4 ) , A41 ) , ( AO (5 ) , A31 ) , ( AD ( 6 ) , A61 ) , 

C( AD( 7) , A71 ) , ( AD( 8) , A81 ) * ( AD (9 ) , A?1 ) , ( AD( 10) , A 101 ) 

CALL DAL ( *0, • C , *0, *0, »C, »0, ♦ 0 , . C, *0, • 0 , *0, ,0) 

CALL SET=>ST ( 4-7 = 001 , .0, A-IDQ0 3, . 0# 4h30C5, • 0, 4MP007, *0j 4HP01 1, *0/ 
C4 u PCl5, .0 ) 

CALL P ESET ( 1000 ) 

C A LL C9 MO UT2 

CALL D G I N I T ( 2 , I V I E w A , 2 0 » I E 9 ) 

CALL 0GIVIK2, IVIE^c/20/ IEP) 

CALL DG I N I T ( 2 , I °L9T , 1 G, IEP) 

CALL CTI\IT( 1, ITEXT,29, IE3 ) 

ENC9CE06, 1; TITLED) 

1 FSP M A T ( ' PEAL T I«E ' ) 

CALL TEXTS ( i , T \ TLEO, 2*, 6 , i , 3 , 3 , I£P) 

E V CSDE(96, 2, TITLED 

2 C5P^AT(' v, 5Tl?Ni ANALYSIS PRSGPa^ ') 

CALL TEXTS ( 1, TITLE 1,24,?, 1,3,3, I ER) 

I = 1000000 
call DELAY 

I V I E'a3 ( 1 ) a I-*EAO( 1 , 10 ) 

IVIEaS(2) =IPACK(0.0, 1 .0*0) 

I V I E v. 5 < 3 ) =IPAC'< (C. 0 , 0 . 5 , 1 ) 

I V I E ^ 3 < 4 ) -IPACX(-1 .0,0. 0,0) 

I V I E a 3 ( 5 ) = I PAC< ( -0 • 5, C • 0 , 1 ) 

I VIE/. 3 (6) ■IP4C , <(0.5»C.0a0) 

IVIEaS(7) =I°aC'<( 1.0, 0*0,1) 

IVIEaP (3 ) =IPAC<(0.75, 1 .0,0 ) 

I V I EX? ( 9 ) = I p ACX ( 0 . 73, 0 .5, 1 ) 

IVIEa5( 10)* !°AC<( 0»5,0. 75,0) 

IV !E*B ( 1 1 ) «I°AC< ( 1.0,0.75, 1 ) 

I V I E a 5 ( 1 2 ) = 0 

CALL GPAPi-a (2, I VIEWS, 12, 1, I ER ) 

IP9 = C 
1 C 9 A = 0 
109 = 0 
\DA =12 
N; A 0 = 1 C 
IN = 0 
X=-100. 

Tir*E=C.O 

A A A A = C . 

P = C.C 

0 = 0.c 
PH9 = 2 • 0 
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DAT£ 



0=32*2 

!-s 1 30 • 

CDX = • 3 
CD Y= 1555 • 0 
LW = 1*0 
L5»1C« 

J= 1 

I Z = 3333 • 

A22 = 1 • Q 
AU= * 45 
AUC9T = .45 
A V = • 0 2 C 
AVD9T= .C20 
AR=0*C100 
A9D° T »C • 0 IOC 
AXO = 25 * 
avQ = 25 * 

AT =3* 75 
AZ-.C05 
IGC = C 
TX = C ,C 
TY = -6000 • 

c5c C a LL DTI'IITdi ITEXT/29 j IE^> 



501 FORMAT ( ' PRESENT ' ,2F10 *2, ' 
C' ' , F6 • 3 ) 


’ , 24 * 


2/ ' 1 


* PS. 2/ ’ 


1 / r 6 • 3 4 


525 F9R^AT( i v-AXIyl^' , 2FlC*2, ' 
C ' , F 6 * 3 ) 


' , F a • 


2/ ' 


/ F 3 • 2 / ' 


1 / r 6 • 3 4 


526 PE-R^AT ( ' v I N I ^ L M 1 / 2 r 10 • 2> ' 
C ' ' , F 6 * 3 ) 


1 ,F4* 


2/ ' ' 


/ ~ 3*2/ 1 


1 t F 6*3/ 


3C7 FORMAT ( 'ORESE vt ' ,"3*2, R1 0*3,' 
CF5.3,' ' , F3 • 3 ) 




1 3 ~ ^ • C 4 


' ' , PS • 3* ' 


f 4 


527 FORMAT ( ' m A X I M ',“3*2,F10*3, ' 
CF5 • 3, ' ' *F5*3) 




1 / r 4 t 2 > 


' ' / p 3 * 3 / ' 


' / 


523 FORMAT ( ' ^ j n I v ij ^ ' , r 8 • 2 » R 1 0 • 3 , 1 
C r 5 • 3, ' 1 , F 5 * 3 ) 




1 / ~ k • 2 j 


' ' , FS • 3/ ' 


* i 


511 F3fi'^A T ( 'PRESENT - , p 7 . 3 , F 1 C • 3 , ' 
C 1 , P6 • 3 ) 


' * p 1 


C • 2 / ' 


' ,2 s.3/ ' 


’ 4 r 6 • 3 m * 


329 F2RN4T( • *AXl'-'J'* ' ,27.3/FlC»3, ' 
C ',26*3) 


' / PI 


C-2, ' 


' / p 3 • 3 / ' 


1 1 r 6 ♦ 3 / 1 


5 3C F?R^aT{ ' v I M M J M ' ,P7. 3/p 10*3/ 1 
0 '/FA. 3) 


'/PI 


0*2/ ' 


' , 2 5*3, ' 


’ 4 F 6 • 3^ f 


5 16 FORMAT ( E6 . 2, ' ',24*2/' 


’ / F4 


*3/ ' 


' ,26*2/ ' 


1 4 “5 • 34 1 


C ' ,27 *2/ ' / ' / F5* 1/ ' ' / F5 . 2, ' 

ENCODE ( 96, 500 , I _ A 1 ) 

500 FORMAT (' 


' / F3 


• 5 / 1 


' / 25 * 2 ) 
=3Sl^I9N ' ) 





CALL TEXie( 1, I LAI, 24, 3, 1/1/3/ IER) 

IP( I23.NE*0)5'JT="jT( 1C1 ) IER, ' = 5S!T;9\" 

ENiC a 0E<96,502, ILA2 ) 



5C2 FOR M A T ( ' X Y 

r°t~ I T U r r /v i ) 

'call T EXTe( 1,?LA2,24, 4, 1,1,3, IER) 
ENCODE <96, 505, ILA3) 

305 FORMAT (' 



ZZ = SI 



VELOC I Tv ' ) 
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DATE 



105 

106 

107 

108 

109 

110 
111 
112 

113 

114 
113 
116 
1 l 7 
112 

115 
120 
121 
122 

123 

124 

125 

126 

127 

128 
125 

1 30 

131 

132 

133 

134 

135 

136 

137 

138 
135 

140 

141 
1 *2 
143 
1 ii4 

145 

146 

147 
143 



20 °T 



CALL TEXTS ( 1, !LA3i24, Uj 1, 1,3* 122) 

I- ( !ER.ME*0)5uT D UT( 1015 122/ ' VELOCITY' 

ENCODE (96/517# ILA4 ) 

517 F82--*AT(' y V w 2 

CP GM 

CALL TEXTOC 1/ ILA4/ 24/ 12/ 1/ 1/ 3/ IE2 ) 

I F ( I ER . NE • 0 ) 0UT D UT (101) 122/ ' _IST vA2 VELOCITY’ 

ENCODE ( 96 / 505/ I L A E ) 

509 FORMAT ( ' ACC2LESA T I SN ' ) 

CALL TrYT6( 1/ ILA5/24/ 19/ 1/ 1/3/ 122) 

I r < IE2.XE.0)9JT°UT( ion IER/ ' ACCEL TITLE' 

ENCODE (96/512/ ILA6 ) 

512 F S 2 M A T ( ' JD9T VD?T LOST 

COT Q C S T ' ) 

CALL TryjSf 1/ ILA6/ 24/ 20* 1/ 1/ 3 / 122 ) 

IF ( I E2 . NE • 0 ) ?‘JT° JT ( 101 ) 122/ 1 ACCEL VA.2' 

ENCODE (96/514/ I L A 7 ) 

514 F 3 p M A T ( ' C0NT29L INPUTS TIME 

Cl GABION DATA ' ) 

CALL TEXTS ( 1/ lLA7/24,27»l/l/3/IE2) 

IF ( I E 2 . \ E ♦ 0 ) S 'J T 3 U T ( 101 ) IER/ ' C S N T 2 3 L HEADING' 

ENCSDE (96/515/ ILA3) 

515 F?R m a t ( ' THRUST 2UCQF2 CELT ELAPSED 

CATE DRIFT 2SLL S^EEO ' ) 

CALL TEXTS ( 1/ILA8/24/2S/ 1/1/3/ 122) 

I F ( I E 2 • N E • 0 ) S L T 3 'J T ( 1CU 122/ ' 

CALL ^ 2 I T E C L S C :< ( 0 ) 

CALL 3TA2TCLSC< 

Id CSNTIN'UE 

CALL ADDA(AO/NAD/AA/NDA) 

CALL 2E A0CL3CX ( N ) 

CALL W2ITECLSC'<(C) 

XI -PS 1/6.28 3 135 

IX = Xi 

AlX=Xi-lX 

IF(AIY.LT.C.C) 3? T S 10 
READ = A I X *360 • 

30 T S 11 
1C A 1X1 = 1.+ A IX 

READ = A 1X1*360* 

11 CONTINUE 

AR ATE =(360* *2 )/6. 233185 
A2SLL=(36C»+ f3L| I )/6.233l85 
AP I T C U = ( 3 60>* THETA ) /6. 28 3 135 



ITCH 



3LS LABEL C3NT2SL' 



14® 






S 3 EED = 


L*< 


>59721 


3 




150 






i f ( m . 


G T - 


.XC )G? 


TO 


700 


151 






HI 


= xc 










152 


700 


IF ( H2 « 


GT, 


. YOGS 


T ^ 


701 


153 






+ 2 


s v r> 










154 


7C 


:i 


IF 


(H3. 


G T < 


-ZZ)GS 


Q 


702 


155 






43 


-11 










156 


7 ( 


:a 


1* 


( 4 4 . 


GT. 


• PS I ) G 


C) T Q 


703 


157 






H4 


II 

u 

O) 

i — a 











N 



FADING 
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DAT 



7C3 


I c (wg.GT. 
45*Pn I 


PHI) GO TO 704 


704 


I r ( H 6 • G T * 
H6»HETA 


theta )3S TO 705 


70S 


I- (H7-GT. 
4 7 = ' J 


J)G9 TS 706 


706 


IFH3.3T. 
1-3 = V 


Vice TO 707 


707 


IF(W9.QT. 
4*9 — A 


W)3e T? 703 


705 


IF (H10.QT 
4 10 = R 


• R ) 39 TS 7C5 


709 


IF ( 1 1 . 3 t 
41 1 :0 


.= )ae to 7io 


710 


I " ( 4 1 2 • 3 T 
412 = 3 


•3)3? TS 711 


711 


I F ( - 1 3 * GT 
4l3=UC5T 


• J 0 S T ) 3 S T2 '’I? 


712 


IF (-14,31 
H 1 4 = V C 5 T 


. V D 5 T ) G 3 T5 713 


713 


IFH15.'3 T 
H 1 5 = ■, 0 9 T 


. W D g T ) 3 9 T 3 714 


714 


l c (H16.3T 
Hl6=RCgT 


•ROST) 39 T9 713 


715 


IF ( u l 7 . RT 
h 1 7=°CgT 


• FOOT ) 39 T9 716 


716 


I F ( h 1 3 • 3 T 
H 1 3 = 0 C 0 T 


•GD° T )3S TS 718 


718 


I r (h3C-.LT 
- 3 C = X 0 


. XO ) 39 T5 71? 


719 


I c (H3! *LT 
43 1 * v 0 


.40)39 T9 720 


72C 


IF ('-32.LT 
H32=ZZ 


. ZZ ) G? T9 721 


721 


IF (’-33.LT 
H 3 3 = p S I 


• PS I ) 3? T9 722 


722 


IF ( 4 3 4 , i_ r 
h3 4 = P>- i 


.pH 1)39 T9 723 


723 


I F ( -35 • LT 
1-35*7 ’-pta 


.TWETAJQg T9 724 


724 


If (436.LT 
436 = 0 


. U ) G5 T g 722 


725 


IFH37.LT 

H37=V 


. V ) G 2 75 726 


726 


IF (433 .LT 
r 3 8 = a 


. w 5 G 9 T9 7 27 


727 


IF ( -39.LT 
439 = 9 


.2)39 Tg 728 


723 


IF ( 44C.LT 
440 sP 


. P ) G? T 9 723 


729 


if<’-41. L t 
44 1=0 


• G ) 39 Tg 730 


73C 


IF ( 4 4 2 . !_ T 


* JD9T ) 39 TS 731 
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DATE 



211 




642=jC9T 






212 


731 


ir(-j45. U T.V93T)G9 


T9 


732 


213 




H 4 3 = 7 C 9 T 






214 


73? 


I r ( 4! 4 4 » 1_ T * W 0 5 T ) 39 


T9 


733 


215 




644= vD9T 






216 


723 


I r (-145. i_T .RD9T ) 30 


T9 


734 


217 




PJ45 = R09T 






213 


734 


IF ( 1J 46«LT.PD5T)38 


T9 


735 


219 




ua6=°D9T 






220 


725 


I tr (H47.LT.009T)G9 


T9 


736 


221 




H47=3D5T 






22? 


736 


C9NIT I \UE 






223 




CELT = ( OFLT + ( \/60 . 


) ) /? 


• 


224 




T I » T I ^E+OELT 






225 


651 


Z=-0.5+a91 






226 




IP ( Z »C-T . .5 ) GP T9 


790 




227 




IF ( Z *LT .-.5)39 T 9 


791 





228 

229 

230 
231. 

232 

233 
23* 

235 

236 

237 

238 

239 
24C 
2^1 

242 

243 

24 A 

245 

246 

247 
2 43 
24? 

250 

251 

252 

253 

254 

255 

256 

257 
253 
259 
26C 
261 
262 
263 



IFCZ.GT. .005)35 T9 792 
IF(Z«LT.-. 005)3? T9 792 
2 = 0.0 
39 T 9 792 

790 Z = * 5 

3? T 9 792 

791 Z= - • 5 

G9 T 3 72? 

792 C9^ T !'IE 

T = A~ *i31 *100. 

I c ( T «3T. 375 • ) 39 T9 793 
!F(T.LT.C.C)G9 T9 79* 

ca tj 705 

793 T = 375 » 

39 T 5 7 Q 5 
7?4 T = 0 ♦ C 

39 T9 795 
795 C5N ,T I - s ;uf5 
998 CONTINUE 

PA AP = 00 . 2C*S I N. { 3 S I ) 

P'~I*F 

T'-h£t a = (j09T 
a-PAA- 
P = p A 
Q=PAA9 
WDoT = P A A P 
P09"=P AAP 
QD6 T = P AA P 
U 3 A 1 1*X*AU 
V=A21*X»AV 
R=A3l*X*AR 

UD6T=( (-l.)*(CDX/18C.)*J*J)4( (T/133. )*C95(Z) )+V*R 
V09T= ( ( T/6 ) *S ! \ ! Z ) ) • ( (CO 7/7! ) *V*A93( V ) ) - j*R 
PD9T=( (COY/IZ) *LW*V*A3S(V) )-( ( J/ I Z ) *!_9 *3 I X ( Z ) )+(a?2*j*V 
SETA s a TA\ ( V/U ) 

X3 = ( ( U*U J ( V *V ) ) **.5 
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La/IZ ) ) 



DAT 



PS I = PS I +DELT * P 

xcdot=u*cos ( =s i ) -v*sr-(Psi ) 

YOD”T=lj + S I N ( P S I ) + V*C9S('S[ ) 

XO=*C+OELT*XCDO T 
YO=vc+OELT*YODD t 
311=UCOT/(A'JDOT*100. ) 

32l=VD0T/ ( AVD0T*100. ) 

331=RC?T/(APD?T*100« ) 

371=X0/( AX0*1C0- ) 

B21 = vO/( AYO^lOO* ) 

IF ( SE V SE SWITCH 1)102,105 
1C2 CALL u 0LD 

CALL ST5PCL OC < 

nFI^F( 6, l00)LF f JD,TTD,TD,T3EA-, J EAD, ATVECDiCCC, TvEC 

100 F 0 R M A T ( 9 F 9 » 4 ) 

CALL ST 4 RTCL0C'< 

CALL COMPUTE 
1 C 5 CONTI v (j E 

IF( 4 10 i • GT * -0 • 5 ) G3 TO 101 

lco= ico+i 

I r ( ICO.E0.1 !G5 TO 653 
I F ( I C G • EC * 3 ) 30 T 3 653 
I c ( ICO. EG. 15)03 TO 653 
I r ( ! CO • EG « 22 ) GO TO 653 
I c ( !C^ .EC. 23 ) 30 TO 653 
I r ( I C Q • E G • 3 6 ) GO TO 653 
GO T 0 107 

6 53 ENCODE ( 96 , 50 1 , ! _ A? ) yC , X 0 , Z Z » =S I , a '-t I , T'NET 4 
6C0 CALL T EXTO( 1 , !l\R,24,6, 1,1,3, IEE) 

IP ! IE D »NE .0 ) OLTP'JT < 101 ) IER, ’ PRESENT' 

107 I F ( IC5.EQ.5JGS TO 108 

q 0 r * 103 

103 ENCODE ( 96,5 25, I_A = ) -2, -1 , *^3, H4, ~<5, 06 
CALL t EXTO (1, I_AP,2^,3, 1,1,3, I E P ) 

109 IF( ICO. EC. 10)30 TO 110 
GO T 0 111 

11C ENCODE ( 96,526, I LAP ) - 31, -'30,-32, -33,63 A, p35 
CALL T EXTO ( 1 , I LAO, 2», 9, 1, 1 , 3, IEE ) 

111 !F< ICO. EC* 2)3? TO 6C1 
IP( ICO. EC. 9)30 TO 601 

I F { ICO.EC* 16)30 TO 601 
[=■( ICO. EQ, 23)30 TO 601 
IP ( ICO. EC. 30) 3^ TO 601 
IF ( ICO .EQ. 37) 30 TO 601 
GO To 112 

6C 1 ENCODE (96,507, ILAP) v ,V, , P , 3 , G 

CALL TEXTO ( 1, I LAP, 2A, 14, 1, 1,3, I 59 ) 

112 IF ( ICO. EC. 15)30 TO 113 
30 TO 1 1 A 

113 ENCODE (96, 327, ILAP ) s-7, h 3, H9#-H C, -ill , ■-'12 
CALL TEXTS ( 1, I LAP, 24/ 16, 1, 1,3, I EP ) 

114 I p ( rCO. EC. 20)30 TO 115 
GO TO 116 
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DATE 



317 


115 


ENCODE (96/ 523, I 


_ A 3 


) 636 > 


637 


, ri 


3? 


,-*39, 


640 


,-.41 


31? 




call texts ( u !lap> 


24,17 


, 1 , 


1, 


3* 


IE?) 






319 


116 


I F < !C a *EQ.3)GS 


TS 


604 














3EC 




I tr ( IC?*EQ.10)GS 


TS 


604 














321 




I r ( ICS. EG. 17)39 


T0 


604 














32? 




I r ( IC2. EC. 24)30 


TS 


604 














323 




IF ( I C ° • E Q » 3 1 ) G S 


T3 


604 














324 




IF ( ICS. EC. 33)03 


T3 


604 














325 




00 T 0 117 


















326 


604 


EXC 5 »CE(96/511i I 


_ A? 


) UDS T 


, VO 


6 T 




D5T , 9 


DST 


, POST, CD0T 


327 




CALL TEXTS! 1/ IL 


Ap, 


24,22 


, 1 , 


ii 


3, 


IEP) 






323 


117 


I* ( IC5.EC. 25)09 


TS 


113 














329 




00 T= 119 


















330 


113 


E'iC a DE(?6>529, I 


LAP ) 613, 


'- J 1 4 


, 6 


15 


,-*16, 


61? 


,413 


331 




call TEXTS* 1, I L 


A P , 


24, 24 


, 1, 


1, 


3 , 


IE? ) 






3 3c 


113 


I F ( I C 2 • E Q . 3 0 ) G 0 


TS 


120 














33? 




GS T 0 121 


















33** 


120 


ENC~0E(56,53C, I L A 3 ) 66?, 


-*43 


, 


4 4 


, ^45, 


61*6, 447 


335 




CALL T EXT0 ( 1, I L A° , 


24,25 


/ 1 , 


1 , 


3, 


IE? ) 






336 


121 


I F ( I C2 . EG ♦ 4 ) 35 


T s 


602 














337 




l c ( ICS. EC* 1 1 ) 00 


TS 


602 














33 5 




IF ( ICS .EG .13)30 


TS 


602 














339 




I r ( IC0.EG. 25)30 


TS 


602 














34C 




I r < IC2. EC. 32)33 


TS 


602 














34 1 




l c ( IC2.EG* 39)03 


TS 


602 














34? 




G0 T 0 923 


















34? 


602 


£MC0DE(96i516i I 


LAP 


)T,Z, 


DEL 


T , 


T I 


'“'E, A? 


ITC 


4 , H £ A o , A 9 A T E , 


344 




r* !"\ 
w S 


















345 




CALL TEX t 5( 1/ ILAP, 


24,30 


, 1, 


1 , 


3 j 


IE?) 






3<*6 


323 


I- ( A 101 «GT. - C *9 


0 ) 3 


S T a 


924 












34- 




u = Q • 


















3^5 




V = C • 



















BETA 



349 



4 = 0 



§ 



353 

351 




355 

353 

35 ^ 



3 = 0 . 

X 0 = 0 . 

v C = C • 



356 

357 

35 3 
35? 
36C 
361 
36c 

363 

364 
36 * 

366 

367 

36 5 



EZ=G. 

= S I =C • 
PHI =0* 
THE T A =0 • 
L'D 5 T = 0 « 
VD S t jO . 

-m 3 5 ” = • 
9 DS T =C . 
PQS^=0. 

cdst=c. 

T 1 7 E = 0 • 

61 = 0 . 

6 ? = 0 - 
63=0 • 



36? 



64 = 0 • 



A P 3 L , S = EE 
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DAT 



q q 



24 



6 03 
58 



n 5 = 0 . 
H6 = C » 
H 7 = 9 • 
h8 = C . 
H9 = C » 



H10 = C 
Hll*Q 
H12 = Q 
H 1 3 = 0 
H 1 4 s 0 
H15*0 
Ml 6 5 0 
HI 7*0 
HI 8 = 0 
h 33*0 
H 3 1 * 0 
H 3 2 = 0 
h 3 3 = 0 
H3 4 = 0 
H35 = 0 
H36 = 0 
H37 = 0 
H3 3 = 0 
m 3 9 * 0 
'-40 = 0 
H4 1 =0 
H.42 = 0 
.-43*0 

H4 4 =0 

- 4 b * 0 
H 4 6 = 0 
H47 s o 



CONTINUE 

I F ( I C9 . EQ « 5 ) Q9 T 9 603 
!F( ICO. EQ. 12 5 39 T9 603 
I F ( IC9. £3*19)39 T9 603 
IF ( !Ce*EQ. 26)39 T9 603 
IF ( IC9. £0.33)39 T 5 603 
IF ( IC9 .EG .40) 39 T9 603 
G9 T9 610 

continue 

continue 



X IX* .45 *C9S ! °h i ) 
Yl V *. 45*3 IN ( PH I ) 
XlXA=(2.5*9)+XlX 
XI 1=2*5 *R 



rlY a=t^ETa-yiv 

X 1 X5 = ( .2 • 5*9 ) - XI X 

Yl Y2=THETA+V1 Y 

XiXC* (2*5*9 ) MO* 1*3 In (?h i ) ) 

y1yc s ThETA+(C»1*C53(F'-*I ' ) 

ICS A* IC9A + l + ( 1 . *!j ) + ( 10. *( IC5A/514 
Yl YD* ( IC9 A/514 , ) *( -1 .0 ) +( THETA ) -0 



) 

4 



55 



DAT 



423 
4 2 A 

425 

426 
4?7 
423 

429 

430 

431 

432 

433 

434 

435 

436 
43 T 
43* 
43? 

44: 

44 1 
44 ? 
443 
44^i 

446 

44 ^ 
443 
445 
453 

45 1 
45 3 
452 

451* 

455 

456 

457 
45 5 
45? 

460 

461 

462 

463 

464 

465 

466 

467 
463 

469 

470 

471 

472 

473 

474 

475 



Y 1 Y E 3 ( I CO A/5 1 4 . ) * ( - 1 .0 + ( THETA ) -0 • 1 4 
xlxD= (X 1X5-0. 05 ) -( ( I CO A/51 4 . ) *0 • 2 ) 
XlXE-X 1X9-0* 05 
XlX^sXlXA+O *05 

XlXG = ( XlXA+C. 05 ) +( ( ICOA/51 4 . ) *0*2 ) 
Y1YP=viya-0« 14 
YlYL sV 1/5-0 * 1 4 

ha ( 0.2/C0S( D HI ) ) 
pU*m*SIN< PH! ) 

YlY2E=vlYE-( I C9 A/2700 • ) 

REARAl=-0.5 
REARA2 = -0 .5 + v 1 Y 
REAPS 1 = -0 » 5 
REA = =2r-0.6-* + YlY 
R E A p c 1 = C * 5 
<EA=C2s-o*5-YlY 
REAPD1=0*5 
PEAP02=-0.64-yiy 
c5 SE3X = 0 * 75+ ( V0/500C . ) 
SESY=0.75*(XC/5000. ) 
!P(Y1YE.5T.-0*6)G3 TO 59 
I C 5 A = 0 

59 CONTINUE 

IV IE.* A ( 1 ) al^EAD (Oi 10 

IV IE.* * ( 2 ) = I PACK (EEAPAi, P£ARA2/0) 

IV IEa a ( 3 ) = IPAC'C ( X1XB/ Yl v 8/ 1 ) 

I V I E a A ( 4) = I c AC<(XlXAyYlYA# 1 ) 

IV I E>. A { 5 ) = IPAC<(REAPC1^9EARC2i 1 ) 

I V I E a A ( 6 ) =I=ACK(REAPAl/ : ?EAPA2 i 1 ) 

IV IE.-, A < 7) = I? AC/ ( PEAP51^9EAR32> 1 ) 

IV IE-vA ( 3 ) = IP4C’<( REARD1/PEARD2j 1 ) 

IV IE- A ( 9 ) = I PACK ( REARC 1> 3 EAPC2> 1 ) 

! r ( = - 3T - .02)35 TO 50 

I- (R.LT.--OP)GS TO 51 

IV IEwa ( 10 ) alaAC* (X1XD/Y1YEE>0 ) 

I V I E > A(ll )*I°AC<(XlXEi/iYEi 1) 

I V I E •. M 1 2 ) = I = A c '< ( X 1 X c > Y 1 Y E > C ) 

I V I E -n A ( 13) = I°AC<(XlXG>YlYEE>l ) 

I V I Ea A ( 1 4 ) » I = AC'< ( 0 . C i 0 • Oi 0 ) 

GS T 0 55 

50 IV IEa'A ( 10) als>AC<(XlXAj YlYAyC) 

I V I E « A ( 11 ) = I p AC< (XlXAi YlYR/ 1 ) 

IV IE. '.A ( 12 ) = I=AC<(REAP01/9EARD2, 1 ) 
IYIEWA< 13) =I=AC<(X1XF. Y1YE/C! 

IV IE.* A ( 14) = IPAC<(X1XG>Y1YEE>1 ) 

GS Te 55 

51 I V I E a A ( 1 0 ) = 1 0 A C < ( X1X5> Y1 y 3iO) 

I V I £/ A ( 1 1 ) = I o A C< ( X 1 X 3 , Y 1 Y-_ , 1 ) 

IV IE*A ( 12 ) = I = ac'< ( RE ARBI/ REAR32 / 1 ) 
IVIEaA( 13) =I°Ac<(X1xE»Y1YEjC) 

IV IE A A ( 14 ) alPAC'< ( XlxDi Y1YEE/ 1 ) 

GO 55 
55 CONTINUE 
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DAT 



IV IE * A ( 15 ) = IP AC* (SESXj SESY, 0) 
IV I F vA ( 16)=IPAC'<(SE£X>SESY/1) 

I V I F A ( 17) = 0 

call graphs ( 2i IVIEWA, I7i2/ IER) 
610 IF ( IC5«LE.40)G9 T? 101 

ice=o 
G9 re ici 

EMD 
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